


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1980 


A signal processing algorithm based on 
multiple microprocessors for an underwater 
acoustic imaging system. 


Vermander, Guy Ronald Arthur 


Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/19053 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist Ser Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ies) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


A SIGNAL PROCESSING ALGORITHM 
BASED ON MULTIPLE MICROPROCESSORS 
FOR AN UNDERWATER ACOUSTIC IMAGING SYSTEM 


Guy Ronald Arthur Vermander 








NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


A Signal Processing Algorithm 
Based on Multiple Microprocessors 
for an Underwater Acoustic Imaging System 


by 


Guy Ronald Arthur Vermander 


December 1980 


Thesis Advisor: G. L. Sackman 
Co-Advisor: U. R. Kodres 





Approved for public release; distribution unlimited 


1197037 1197036 





SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 


REPORT DOCUMENTATION PAGE 


2. GOVT ACCESSION NO 





READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


3. RECIPIENT'S CATALOG NUMBER 













5S. TYPE OF REPORT & PERIOD COVERED 
Master's Thesis 
December 1980 


6. PERFORMING ORG. REPORT NUMBER 


4. TITLE (end Subtitie) 











A Signal Processing Algorithm Based on 

Multiple Microprocessors for an Underwater 
Acoustic Imaging System 
- AUTHOR(e) 






8. CONTRACT OR GRANT NUMBER(e) 







Guy Ronald Arthur Vermander 






10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 


12. REPORT DATE 
December 1980 
13. NUMBER OF PAGES 


85 


15. SECURITY CLASS. (of thie report) 


Unclassified 


- PERFORMING ORGANIZATION NAME ANO ADORESS 
Naval Postgraduate School 
Monterey, California 93940 

















- CONTROLLING OFFICE NAME ANDO ADORES3 
Naval Postgraduate School 
Monterey, California 93940 
















- MONITORING AGENCY NAME @ ADORESS(!f different fram Controlling Office) 


Naval Postgraduate School 
Monterey, California 93940 








Se. DECLASSIFICATION 


( OOWNGRAOING 
SCHEDULE 


. DISTRIBUTION STATEMENT (of thie Report) 


Approved for public release; distribution unlimited 





17. OISTRIBUTION STATEMENT (ol the ebetract entered in Block 20, ti dillferent from Report) 


- SUPPLEMENTARY NOTES 









- KEY WORDS (Continue on reveree aide tf neceeeary and tdentify ty bieck number) 


Signal Processing, Microprocessors, Underwater Acoustics, Acoustic Imaging 





20. ABSTRACT (Continue on reveree side if neceeeary and identify by bleck mumter) 


An algorithm has been designed to provide near real-time location and 
classification information to an operator of an underwater acoustic imaging 
System. The overall system to detect objects buried up to five meters in 
unconsolidated marine sediments consists of a specialized projector co-located 
with a five meter receiving line array. This system is described in basic 
terms together with the principles of the signal processing that extracts 
information from the scattered acoustic signals. Using simulated objects, the 







DD orn, 1473 cortion oF 1 Nov 68 1s OBSOLETE 


JAN 73 
Ss N = — 01 eee ee ee eee 
_— ' SECURITY CLASSIFICATION OF THIS PAGE (When Dete Entored) 


1 








ee 
smeTY CL EDIFPIC ATION OF 7m1S PAGE(Wren Nore Entered. 






angular resolution and imaging performance of a number of system designs 
were evaluated. To obtain near real-time execution, the algorithm was 
optimized and partitioned using parallel, pipeline, and double buffering 
techniques for independent but synchronized operation on multiple micro- 
processors. The feasibility of the design approach was experimentally 
emonstrated using four single board computers in a microcomputer develop- 
nent system. 


OQ. 


For 1473 
| Jan ; 2 
0102-014-6601 SECUMITY CLABSIFICATION OF THIS PAGECHROn Date Enioved) 





Approved for public release; distribution unlimited 


A Signal Processing Algorithm 
Based on Multiple Microprocessors 
for an Underwater Acoustic Imaging System 


by 
Guy Ronald Arthur Vermander 
Captain, Canadian Forces 
B.S., Royal Military College of Canada, 1966 


Submitted in partial fulfillment of the 
requirements for the degrees of 


MASTER OF SCIENCE IN ENGINEERING ACOUSTICS 
and 
MASTER OF SCIENCE IN COMPUTER SCIENCE 
from the 


NAVAL POSTGRADUATE SCHOOL 
December, 1980 





ABSTRACT 


An algorithm has been designed to provide near real-time location 
and classification information to an operator of an underwater acoustic 
imaging system. The overall system to detect objects buried up to 
five meters in unconsolidated marine sediments consists of a spe- 
cialized projector co-located with a five meter receiving line array. 
This system is described in basic terms together with the principles of 
the signal processing that extracts information from the scattered 
acoustic signals. Using simulated objects, the angular resolution and 
imaging performance of a number of system designs were evaluated. To 
obtain near real-time execution, the algorithm was optimized and parti- 
tioned using parallel, pipeline, and double buffering techniques for 
independent but synchronized operation on multiple microprocessors. 

The feasibility of the design approach was experimentally demonstrated 


using four single board computers in a microcomputer development system. 
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I. INTRODUCTION 


The specific situation which led to this research is the desire 
to detect, locate, and classify objects buried in unconsolidated 
marine sediments up to a depth of 5 meters. To accomplish this an 
acoustic imaging system is needed which can be mounted on an under- 
water vehicle maneuvering close to the ocean floor. The problem is to 
penetrate the soft sediment acoustically, detect the scattered signals, 
and derive from them information as to the nature of the scattering 
objects. -With sufficient range and angular resolution it would be 
possible to adequately localize and classify an object. Furthermore, 
it is necessary to obtain this information as rapidly as possible: 
ideally a display would provide a "real-time" representation of buried 
objects to an operator as a section of sediment is being searched 
acoustically. 

Although a number of seetnigues of acoustic imaging exist [1], 
a system similar to a high resolution sonar has been proposed which 
would have high resolution using very short baseband pulses and which 
would operate with dynamic focussing deep in the near field of the 
receiving aperture. This system would involve a specialized para- 
metric projector co-located with a receiving horizontal line array of 
hydrophones. It would be mounted on an underwater vehicle operating at 
about 5 meters above the ocean floor, and would provide acoustic 
penetration of the marine sediments to depths of about 5 meters. The 
objective of this thesis was to design the signal processing algorithm 


for such a system. 





In the course of the development of the algorithm the effect of a 
number of system parameters on angular resolution and imaging perform- 
ance was evaluated for a number of possible implementations of the 
acoustic system. This led to major decisions as to the form of the 
algorithm and the possible acoustic imaging system. 

Since an acoustic imaging system mounted on a maneuvering under- 
water vehicle would be required to provide an operator with a near 
real-time display of buried objects, the algorithm was redesigned for 
increased execution efficiency and implementation on multiple micro- 
processors. Optimization of the algorithm by removing time consuming 
operations, together with incorporation of parallel and pipeline 
processes and double buffering concepts were the approaches taken. 
Experiments were aeraucted using a number of single board computers to 
test and demonstrate the feasibility of near real-time execution of the 
algorithm. 

Section II of the thesis describes the acoustic imaging system and 
the basis for the signal processing algorithm. The details of the 
algorithm and the angular resolution and imaging performance evaluations 
are explained in Section III. The design of the algorithm for near 
real-time execution on multiple microprocessors and the experimental 
demonstration of feaSibility are contained in Section IV. The resulting 
conclusions are then given in Section V, followed in Section VI by 
recommendations for future work that could produce a special purpose 


processor for the acoustic imaging system. 





II. THE PROPOSED ACOUSTIC IMAGING SYSTEM 


A. INTRODUCTION 

An acoustical imaging system mounted on an underwater vehicle wil] 
consist essentially of three major parts: a sophisticated projector, a 
receiving array, and a signal processor and display system. The thrust 
of this thesis has been to develop the signal processing algorithm. To 
understand the underlying concepts, the potential form of the projector 
and receiving array is briefly described, followed by the proposed 


algorithm. 


B. TRANSMISSION 

Acoustic waves: can penetrate sediment to useful depths with 
acceptable absorption loss if the frequency is low enough; less than 
about 20 kHz. Very short acoustic pulses (50 to 100 microseconds) are 
required for adequate range and angular resolution. For the concept 
that has been developed in this thesis it has become apparent that it 
would be necessary to have effectively a point source for the projector. 
Baseband pulse generation using a parametric source [2] appears to be a 
promising technique for this application. Such a projector would form 
a pencil beam in the collimated region of its near field that could be 


Scanned to produce a tomographic image of a region in the sediment. 


C. RECEPTION 
1. Physical Configuration 
The scattered acoustical energy would be received by a horizon- 


tal line array (mounted on the underwater vehicle) of a large number of 
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equally spaced hydrophones. The projector would be located at the center 
of the line array. The length of the array would be a compromise between 
the need for a large aperture (on the order of 10 meters for angular 
resolution) and the physical reality of deploying and maneuvering an 
underwater vehicle with the array attached (5 meters probably the 
maximum). 

Moving over the area of interest the underwater vehicle would 
maintain an altitude of approximately 5 meters above the sediment. As 
the desired depth of penetration into the sediment is also 5 meters the 
total range of interest lies between 5 and 10 meters from the center of 
the array. A 60 degree sector (between array bearings of -30 and +30 
degrees from the normal to the line array) provides horizontal 
coverage of 10 ferers in the sediment at 10 meters range from the 
center of the line array (Figure 1). 

2. Wavefront Curvature 

For ranges between 5 and 10 meters, operation of the acoustical 
imaging system is deep in the near field of the line array. The far 
field begins at a range of D¢7) where D is the diameter of the aperture 
and \ is the wavelength [1]. Using a 50 microsecond pulse and a 5 
meter line array this is on the order of 300 meters. The near field 
operation is evident by the marked curvature of the wavefronts of the 
acoustic energy as they arrive at the line array after scattering from 
points within the area of interest. Since the curvature depends on the 
range to a scattering point an imaging processor requires some form of 


dynamic focussing. 
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Figure li. The physical relationship between the 
five meter line array of hydrophones 
and the sector of interest. — 
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Figure 2 depicts four such wavefronts (that originated at the 
indicated bearings and ranges) as a function of time of arrival and 
position along a 5 meter line array. Using 1500 meters per second as 
the soundspeed, the time of the earliest arrival of interest is 6.220 
milliseconds after pulse transmission, and the latest is 14.304 milli- 
seconds. 

3. The "Trace" of a Point Scatterer 

The line array consists of a large number of equally spaced 
discrete hydrophones (each considered as a point) which detect the 
pressure amplitude of an arriving wavefront. The time sampling of the 
outputs of 16 hydrophones in a 5 meter line array resulted in a quanti- 
zation of the wavefront arrival times as shown in Figure 3 (in which 
an acoustic pulse width of 50 microseconds and sampling every 50 micro- 
seconds, beginning 6 milliseconds after pulse transmission, was used). 
For example, considering the wavefront originating at a scattering 
point located at a bearing of zero degrees and a range of five meters, 
the hydrophone sampling produces the following time indices for the 
wavefront arrivals: 21, 19, 17, 16, 15, 14, 13, 13, 13, 13, 14, 15, 
fone, 19, 21. 

The sequence of time indices representing the arrival of a 
wavefront, n(r,@) = (nj Mo.. np) for K sensors, is unique for the 
scattering element at a range r and bearing @: this sequence is termed 
the "trace" of that particular point (related to the concept of the 
“Time Delay Trace Function" [3] ). The traces for all points at the 
range and bearing increments in the area of interest can be pre-calcula- 


ted and stored for later retrieval. 
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Figure 2. The arrival times of four wavefronts as 
a function of position along a five meter 
array. 
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D. SIGNAL PROCESSING 
1. The Basic Concept 

In developing the algorithm two situations were examined: 

a. Insonification of the entire sector by a short pulse 
(broad beam in the plane of the line array but very narrow beam in the 
plane perpendicular to the line array) followed by processing returns 
from the entire sector, and 

b. Transmission of the acoustic pulse in a narrow beam and 
Sweeping across the sector with signal processing occurring at each 
bearing increment. 

Essentially, after a suitable gate time after pulse trans- 
mission, each sampled hydrophone output provides a time series record 
of the scattered acoustic amplitudes arriving at the line array. 
Storing the results (after A/D conversion) at each sample time, n, of 
all K hydrophone outputs causes the memory to contain a time represen- 
tation of the scattered acoustic pressure field, D(n.k), k=1..K, n=l..N 
where N is the maximum sample time. The presence of a scattering point 
will be indicated by the presence in memory of signal amplitudes 
occurring at the appropriate times in each hydrophone record (a pattern 
representing the wavefront curvature). 

The sequence of time indices which uniquely represent a scat- 
tering point are available as the pre-calculated trace for that point. 
An estimator of the presence of a pattern is the summation of the signal 
amplitudes at the appropriate time indices or addresses in each of the 


hydrophone records. For example, if the signal amplitude at the 
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wavefront (from a point at r,@) was considered unity, then the total 


Summation along n(r,@) would equal K, the total number of sensors: 


K 
A(r,0) = >. D(ny 4k) = K 
KI 


where n, are elements of n(r,@). 


As a result of noise and errors in geometry the sum would be less than 
K but would nevertheless have a maximum 1f a scattering point was 
present. 

In sunmary, what is required is to form in memory a represen- 
tation of the scattered acoustic field, search through it systematically 
for wavefront patterns by summation along pre-calculated traces for the 
ranges and bearings of interest, and recording the resultant total 
amplitudes as a function of range and bearing. A decision process 
would then follow to determine which amplitudes represented point 
scatterers, as opposed to solely noise, and then these could be dis- 
played forming the acoustical image of the insonified region in terms 
of range and bearing or rectangular coordinates. 

2. Data Structures 

To visualize the signal processing operation it is useful to 
discuss it in terms of the data structures involved: the “Indata" 
array, the “Trace” matrix, and the "Amplitude" array as shown in 
Figure 4, 

The Indata array contains the time record of the scattered 


acoustic field. It can be visualized as two dimensional in memory: 
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Figure 4, The basic data structures 


17 





i.e. D(n,k) where n is the time index, and k is the hydrophone number. 
Each row would be formed by the simultaneous sampling of all hydrophones 
at intervals equal to the pulse width of 50 microseconds, beginning 
after a gate time of 6 milliseconds, and storing the results. This 
would continue for 167 samples until the last possible arrival time of 
a signal of interest (time indices 0 to 166). 

The Trace matrix, T(@,r,k) contains the traces n(r,@) for 
points in the area of interest as identified by the intersections of all 
range and bearing increments. It can be visualized as three dimen- 
Sional where the first dimension is that of bearing, the second is that 
of range, and the third is the hydrophone number. For a particular 
range and bearing, the elements in the third dimension contain the time 
index sequence n(r,@) that is the pre-calculated trace: i.e. n, = 
ieosr,k). 

The dimensions for the Amplitude array are range and bearing. 
Each element, A(r,8), is Ne total sum of the contents of those 


elements from the Indata array identified by the trace n(r,6): 


k 
i.e. A(r,@) =)" D(n,,k) 


k=] 


where ny, are elements of n(r,8) 


Scattering points in the acoustically scanned region will be evident 
by larger amplitude elements in A(r,6) and the distribution of these 


will comprise an image of any objects present. 
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III. THE DEVELOPMENTAL ALGORITHM 


A. INTRODUCTION 

The signal processing algorithm was initially developed using a 
Hewlett-Packard 9845B desktop minicomputer (Figure 5). This machine 
uses HP enhanced BASIC [4} and was equipped with 192K of storage, two 
peripheral flexible disk drives (HP9885M and HP9885S) , a Graphics ROM 
(98437B), and a CRT Graphics Memory (98470B). This proved to be a most 
convenient tool for the development of the algorithm, evaluation of 
the influence of the major parameters involved, and for the graphical 
display of simulated objects, all of which led to major conclusions and 
decisions as to the form of the algorithm and the possible physical 


acoustic imaging system. 


B. DESCRIPTION OF THE ALGORITHM 
1. Parameters 

The parameters involved in the algorithm (Appendix A) were 
those which describe the physical implementation of an underwater 
acoustic imaging system and the area or sector of interest within the 
sediment that is to be insonified. These parameters were varied to 
evaluate their influence on the performance of the acoustical imaging 
system. 

The lengths of the horizontal line array used in the develop- 
mental algorithm were 10 meters and 5 meters; indicative of the desire 


for the largest aperture possible for angular resolution, and the most 
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probable realistic maximum for a physical structure to be attached to 
a maneuvering underwater vehicle. Equally spaced hydrophones along 
the line array numbering 16 and 32 were used, as were pulse widths of 
100 and 50 microseconds. The sound speed was taken as 1500 meters 
per second throughout the transmission path of the acoustic pulse. 

The sector of interest within the sediment was defined as 
lying between the ranges of 5 and 10 meters from the center of the line 
array and between bearings of -30 and +30 degrees from the normal to the 
line array. Range increments of 0.1 meters, and bearing increments of 
two degrees were chosen. This corresponded to 51 range indices from 0 
to 50, and 31 bearing indices from -15 to 15. Bearing increments of 
one degree were used in the evaluation of array bearing resolution. 

2. Calculation of Traces 

The trace, n(r,@), for a scattering point at r,6 consists of 
the time indices representing the arrival of the scattered wavefront 
at each of the sensors in the line array. With the projector at the 
center of the line array, the distance from the projector to each 
sensor was calculated. A simple trigonometric relation was used to 
calculate the travel time for an acoustic pulse from the projector to 
a point P({r,@) and then to the k'th sensor. 

The geometry of two cases 1s illustrated in Figure 6 and 
Figure 7. The first case, Figure 6, applied when the k'th sensor was 
past the array midpoint; i.e. k > K/2 where K is the number of sensors 
in the line array. For this situation the trigonometrical relation 
providing the distance d, from the k'th sensor to the point P(r,@) in 
terms of the range r, bearing 8, and sensor to projector distance b,, is 


as follows: 
21 





] k K sensors 


Tiné array 





P(r ,@) 


Figure 6. The geometry for k > K/2. 
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Figure 7. The geometry for k < K/2. 
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Za 2 ‘Aiea a 
dy by +r 2 bY r cos(90-8) 


Similarly, for the case k < K/2 (Figure 7), the distance d, is 


obtained from 


d, @ = Pea + re - 2 by r cos(90+8) 


The total travel distance for a scattered pulse, r + d,, divided by 
the sound speed, gives the time that the leading edge of the pulse 
arrives at the k'th sensor after transmission. The period of 6 milli- 
seconds is subtracted from the travel time, which is equivalent to 
Starting data acquisition after gating out the first 4.5 meters in 
range from the center of the line array. 

The integer portion of the quotient of this time and the 
sampling interval (equivalent to the pulse width) provides the time 
index ny. Computed for each of the K sensors this results in the 
sequence n(r,@) = {n,} where k = 1..K. This computation was repeated 
at all ranges for each of the positive bearings. Symmetry about the 
normal to the line array permits the calculations of travel times for 
positive bearings to suffice. Negative bearings use the mirror image 
of the trace sequence; i.e. n(r,-9) = {n,} where k = K..1. 

3. Simulation of Objects 

In order to develop the algorithm and in the absence of real 
hydrophone inputs, sampling, and I/0 operations it was necessary to 
Simulate the Indata array D(n,k). An idealized situation consisting of 


zero noise and unity signal amplitude at the wavefront of a scattered 
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pulse was considered. In other words no attempt was made to model a 
noise environment or the reflectivity of a buried object. 

A point scatterer was simulated by using the pre-calculated 
trace, as stored in the Trace matrix T(8,r,k), to identify the time 
indices n of D(n,k) that represented the wavefront arrival at each 
sensor. Unity signal amplitude was inserted into these elements. For 
the narrow acoustic beam case only a Single point scatterer was assumed 
to exist at each range increment along the bearing of the beam. It was’ 
more complicated in the case where the entire sector was assumed insoni- 
fied. In this case the acoustic field would be composed of the inter- 
ference pattern of numerous wavefronts if the sector contained a 
distributed object composed of a number of point scatterers. 

For the latter case two input situations from the sensors were 
Simulated: a nonlinear input, and a linear input. The nonlinear 
Situation would arise if the outputs from the hydrophones were ampli- 
fied, passed through a thresholding comparator circuit, and a binary 
input presented to the computer I/0. The unity amplitude mentioned 
above simulated this situation. If the sampled outputs from the 
sensors were subjected to A/D conversion prior to storage in memory, 
then a more accurate or linear representation of the actual pressure 
field would be available. This was simulated by the summation of the 
wavefront amplitudes in the appropriate elements of D(n,k) as each of 
the points comprising a distributed object was considered. 

The angular resolution of the imaging system was evaluated by 
Simulating discrete points in the area of interest. Imaging perform- 


ance was studied using simulated objects composed of a number of 
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scattering points. The resulting Indata arrays were then operated 
upon by the algorithm to form the Amplitude array which represented the 
acoustic image of the insonified sector. 

4. Forming the Amplitude Array 

The magnitude of each element in the Amplitude array, A(r,@), 
is an estimator of the presence of a point scatterer at the range r, 
and the bearing 8. Each element results from the sum of those elements 
in the Indata array, D(n,k), identified by the trace n(r,9) stored in 
the Trace matrix, T(@,r,k). 

For the case of the entire sector being insonified, the Ampli- 
tude array was formed by summing along the traces for all ranges at 
each bearing in turn for the entire sector. For the narrow beam case, 
only the traces corresponding to all ranges at the beam bearing were 
utilized, followed by the acquisition of a new Indata array for the next 
bearing. Then the traces for all ranges at that bearing were used to 
guide the summation process, continuing until the narrow beam had 
Scanned the entire sector and the contents of the Amplitude array 
represented the presence of any point scatterers within the sector. 

5. Graphics Display 

To evaluate angular resolution and imaging performance of the 
proposed system the graphics capability of the HP9845B was utilized. 

It was possible to draw graphs representing the received amplitudes 
from point scatterers and thereby obtain the system's angular resolu- 
tion. Both a polar or range-bearing grid and a rectangular grid were 


devised to provide a means of referencing a displayed image. It was 
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decided that the rectangular grid would be more practical and 
meaningful to an operator of an acoustic imaging system. 

Due to the choice of 0.1] meter increments in range and two 
degree increments in bearing, a single point scatterer was represented 
by graphically "filling-in" an area on the screen equivalent to these 
dimensions and centered at the range and bearing of the point. This 
was termed a "target" in the developmental algorithm. A threshold 
was used to determine when the magnitude of an element in the Ampli- 
tude array indicated a target present. An object was simulated by 
uSing a number of these discrete points which the graphics display 
then presented as a number of adjacent targets. This gave the image 
outline a stepped appearance rather than the smooth lines of a real 


object. 


CREO ULTS 
1. Range Resolution 

The use of an extremely short pulse of 50 microseconds has 
inherent range resolution much greater than required or utilized in the 
development of the algorithm. It was decided that a range resolution 
of 0.1 meters was adequate to localize and classify large buried 
objects. Accordingly, it was the range increment chosen for the pre- 
calculation of all the traces, and the subsequent processing and 
display of the acoustic image. An additional consideration was the 


memory size available to store the Trace matrix. 
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2. Bearing Resolution 


For the case where the entire sector is insonified by one pulse, 
the effects of a number of parameter changes on the angular resolution 
of the acoustic imaging system were obtained. The bearing resolution 
was determined by the ability of the system to distinguish between two 
discrete points at the same range but at slightly different bearings. 
The two points were simulated as being in the vicinity of a bearing of 
zero degrees, and then a bearing of thirty degrees, and at ranges of 
five meters and ten meters. 

The algorithm processed the resulting Indata array in each 
case and the bearing resolution was determined from the resulting 
amplitude distribution in bearing at the range of the points. Figure 
8 gives an example of the amplitude distribution for two points at a 
number of bearing separation angles. The amplitude distributions were 
Superimposed as the bearing separation between the points was increased 
from three degrees to six degrees in one degree increments. The separa- 
tion in bearing that resulted in 3db and 6db dips between the peaks in 
the amplitude distribution were recorded and are shown in Table 1. 

This was done for both the nonlinear input case and the linear 
input case, and for sixteen and thirty-two hydrophones in the line 
array. It was observed that little difference in bearing resolution 
resulted as could be expected for the idealized zero noise situation 


that was being simulated. 


27 








"ACAAR DULL] 
Aaya ¢ We wory abues suaqzau OQ 2e sarbue uojzeuedas Butueag 
$O saquinu e 72 SzULOd OM 4OJ SUOLYNGLUYSEP apnqytdwe auy 


(2350) ONISWSE 


6 8 é 3 S bP 





°g aunbry 


=, 
— 


SONLT dW LiNdNI 20 WIS 


28 





Sdaqaweued wWazZsds 
JO UOLZIUNJ & SB UOLYNL[OSaYy BHuLueag “| JISVL 








(3ndut 
AeauL|[) 

Aesse 
Jaya ¢ 








(qndut 
APaut[uoU) 

Aeuse 
Jaya ¢ 








(andut 
APoutL[LUoU) 

Aesse 
Jazow OL 











(saaubap) 
butueag 


(suaqzaow) 
abuey 


( spuodasoud tw) 
UIPEM 9S5[Nd 





(saaubap) uozqnposay 





29 





3. Simulated Images 


Three large objects, each considered as being composed of a num- 
ber of point scatterers, were simulated in the sector or the field of 
view of the system. One was small and rounded as if it were a cross 
section of a long cylindrical object. The other two were long narrow 
objects; one was in a horizontal position, and the other at an inclined 
aspect angle. Figure 9 shows the positions and orientations of the 
Simulated objects. 

For the case of a pulse insonifying the entire sector and 
with nonlinear inputs from the receiving hydrophones, the maximum 
amplitude for any element in the Amplitude array was equal to the number 
of hydrophones. The threshold for the display of the resulting targets 
was varied to see the spatial distribution (in range and bearing) of 
the amplitudes of the image elements. Figure 10 shows the smearing of 
the objects resulting from the lack of angular resolutions; i.e. the 
skirts of the amplitude distribution in bearing of each scattering 
point were being detected. 

Similarly Figure 1] shows the results of having linear inputs 
from the hydrophones. The largest amplitude elements were considerably 
larger than the nonlinear case, and the contours of the amplitude dis- 
tribution in the image were much steeper. It was observed that the 
size, shape and aSpect angle of the objects had a substantial effect 
on the resulting amplitude distribution in the image. This would be 
beneficial for the identification or classification function of an 
acoustic imaging system particularly if it provided an operator con- 
trolled threshold for the display or if the intensity of the display 
depicted the amplitudes of the image elements. 
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For the case of a narrow two degree beam width system that 
transmits a pulse at each bearing increment, the resulting image after 
the sector was completely scanned is shown in Figure 12. It was ob- 
served that (as expected) the image exactly duplicated the simulated 


objects. 


D. CONCLUSIONS 

The developmental algorithm was used to study some of the effects 
of the system parameters on bearing resolution, and the imaging per- 
formance of possible implementations of an acoustic system that used 
the proposed signal processing algorithm. It was concluded that the 
acoustic imaging system should be a narrow beam system that scans the 
sector by transmitting a very short pulse at each bearing increment, 
provides linear inputs to memory from the hydrophone line array, and 
provides an operator controlled threshold for the image display. Al- 
though the benefits of the latter two were not evident in the simulated 
idealized narrow beam images, it is felt that they would be significant 
in a real system where noise and objects with different scattering or 
reflective properties were present. This is also the basis for 
requiring a line array as a large receiving aperture, together with the 
proposed signal processing algorithm, since objects with specular 
reflective properties would not necessarily scatter the acoustic energy 


in the direction of the projector. 
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TV. THE NEAR REAL-TIME ALGORITHM 


A. INTRODUCTION 

An acoustic imaging system mounted on a maneuvering underwater 
vehicle would be required to provide an operator with a near real- 
time display of any objects buried in the sediment as the vehicle 
conducts a search pattern over the ocean floor. Although the signal 
processing algorithm that was developed was basically simple computa- 
tionally, it was nevertheless very time consuming and would have been 
inadequate in an operational environment. Consequently, it was re- 
designed for increased execution efficiency and implementation on 
multiple microprocessors while maintaining the basic signal processing 
approach that has been described. Optimization of the algorithm by 
removing time consuming operations, together with incorporation of 
parallel and pipeline processes and double buffering concepts were the 
approaches taken. Experiments were conducted using a number of single 
board computers to test and demonstrate the feasibility of the rede- 


Signed algorithm for near real-time execution. 


B. THE TIME PROBLEM 

An underwater vehicle acoustically scanning the sediment and moving 
at 1 meter per second would probably require a complete sector scan at 
least once per second as a minimum adequate display rate; i.e. a two 
degree beam width, approximately 0.34 meters wide at 10 meters range, 


would scan a 60 degree sector perpendicular to the direction of travel 
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once every meter of forward travel. At each bearing increment data 
acquisition would begin 6 milliseconds after pulse transmission and 
continue for 8.5 milliseconds. Thus the processing of the Indata 
array can begin 14.5 milliseconds after pulse transmission, and would 
need to be completed in about 15 milliseconds for each of the bearing 
increments if a scan were to be completed in about one second. It was 
obvious that the developmental algorithm would require considerable 


changes for an efficient microprocessor implementation. 


C. ALGORITHM OPTIMIZATION 

The algorithm was first examined with the objective of removing 
time consuming operations such as multiplication. These occurred in the 
addressing of the multi-dimensional arrays in the algorithm which were 
stored in row At Onder Four additional one dimensional arrays were 
introduced which provided a table lookup of the row addressing. For 
example, the addresses for the bearing dimension of the Trace matrix, 
T(8,r,k), were stored in B_array(b), and those for the range dimension 
in R_array(r). Instead of the addressing for a single element of 
T(6,r,k) requiring three additions and two multiplications, it was 
addressed using T(B_array(b) + R_array(r) + k); i.e. five addition 
operations. This was done for each of the multi-dimensional arrays, the 
Indata array D(n,k), the Trace matrix T(9,r,k), and the Amplitude array 
A(r,@), thus reducing the data structures in the algorithm to strictly 
one dimensional arrays. 

With the intention of testing the feasibility of a microprocessor 


implementation of the algorithm using a number of Intel SBC 80/20 Single 


I) 





Board Computers (Intel 8080 CPU) [5], the algorithm (Appendix B) was 
written using the high level language PL/I-80 [6]. Executing the 
program on one SBC 80/20 it was found that after the pre-calculation of 
the traces in T(9,r,k) the signal processing aspect of the algorithm 
required 6.45 seconds to complete one sector scan (approximately 210 
milliseconds at each bearing, and not including the time required for 


pulse propagation and data acquisition). 


D. ALGORITHM PARTITIONING 

The structure of the signal processing algorithm was examined in 
order to identify functions that could be accomplished by parallel and 
pipelined processes executing on a number of microprocessors for in- 
creased time efficiency. The concept of double buffering was also 
incorporated to reduce execution time. 

1. Parallel Processing 

The processing carried out at each of the range increments on 

the Indata array, for the bearing of the transmitted pulse, is indepen- 
dent of any of the other ranges. This is the summation that results 
in an element of the Amplitude array A(r,9). It was decided to parti- 
tion this function into separate processes that were identical except 
for the scope of the range indices. These were called the Sum_ampli- 
tude processes and were executed in parallel on a number of micropro- 
cessors. Using three microprocessors for example, process 34 per formed 
summations along the traces for the range indices 0 to 16, process S 


Z 
used range indices 17 to 33, and process 53 covered range indices 34 to 
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50. In this manner all ranges were processed in one third the time 
that was required by a single process executing all range increments. 
2. Pipeline Processing 

The function of data acquisition is called the Indata process 
and was seen to be amenable to pipeline processing with the Sum_ampli- 
tude processes. This process could control the timing of pulse trans- 
mission, sampling the hydrophone outputs, and storing the results into 
two buffer arrays. The Indata array of the developmental algorithm 
became two arrays called the A_indata array and the 8 indata array. 
After the A_indata array was filled with the data from a particular 
bearing, the Sum_amplitude processes were permitted to begin their 
operations on the A_indata array. Simultaneously, the Indata process 
proceeded to fill the B indata array with the data from the next 
bearing increment. When the Sum amplitude processes had finished with 
the A_indata array, they processed the B_indata array. The Indata 
process could then utilize the A_indata array for the next bearing and 
SO On in a pipeline fashion. Figure 13 illustrates the relationship 
between the processes. 

3. Semaphores 

Communication between the processes was accomplished with 
binary semaphores. The Indata process set the semaphores A_data and 
B_data whenever it was utilizing either the A_indata or B indata array. 
Similarly, each of the parallel Sum_amplitude processes set a semaphore 
depending on the data array it was processing. Thus for three parallel 


Sum_amplitude processes six semaphores were utilized: A_ampl, A_amp2, 
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A_amp3, and B_amp], B_amp2, B_amp3. The Indata process waited for each 
of the Sum_amplitude semaphores for one of the data buffer arrays to be 
unlocked before it could. proceed to fill the buffer with new sampled 
data from the hydrophone outputs. Similarly, each of the Sum ampli- 
tude processes waited on the A _data or B_ data semaphore before it 

began processing the data. In this way completely independent but 
synchronized operation of all the processes was achieved. The opera- 
tional form of the algorithm in Appendix C illustrates the parti- 


tioning of the processes and the use of the semaphores. 


E. EXPERIMENTAL DEMONSTRATION 
1. Equipment 
To test and demonstrate the feasibility of near real-time 
execution of the redesigned algorithm an Intellec 800 Microcomputer 
Development System (MDS) [7], shown in Figures 14 and 15,with the 
CP/M monitor control program [8] was used. The MDS is capable of 
being used as a partial system simulator to check out software execu- 
ting from RAM in the Intellec MDS. A total of 64 kilobytes of memory 
was available and four Intel SBC 80/20 Single Board Computers [5] were 
inserted into the MDS motherboard as shown in Figures 16 and 17. 
2. The Processes 
The contents of the Trace matrix, the pre-calculated traces, 
were computed using the program of Appendix B, and stored in a floppy 
disk data file for subsequent loading into memory. The redesigned 
“operational” algorithm (Appendix C) was written as if it were to be 


executed on a Single processor, and consisted essentially of a Setup 
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The system used for the experimental demonstration 


Figure 14. 
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procedure that calculated the contents of the address lookup arrays, the 
Indata procedure, and the Sum_amplitude procedure. This was compiled and 
linked so the program started at location 3000 Hex. The choice of 3000 H 
as the initial address for the program was based on the fact that addres- 
ses 3000 H to 3EFF H were available as private RAM on each of the SBCs. 
All the data structures and variables such as the semaphores that were 

to be global to all the processes were located in common memory above 
4000 H that was accessible by all the SBCs. An empty dummy array of the 
necessary size was declared to ensure the assignment of common memory 
addresses to the global variables. 

Using the MDS and.the Dynamic Debugging Tool (DDT) program [9] 
the program was altered at the machine code level (hex representation) 
by removing the appropriate call statements and creating two versions: 
one which executed only the Indata procedure, and the other the Sum_am- 
plitude procedure. In this way the two processes stored in the private 
memory of Separate SBCs could operate independently but in a syn- 
chronized manner using the semaphores in common memory. The data that 
was operated on, as well aS all the data structures, were available in 
common memory. 

The Sum_amplitude version was then amended using DDT so that all 
the local variables in the procedure, which the compiler had assigned to 
memory locations in the common area, were readdressed so they were 
within an SBC's private memory. The locations assigned to the dummy 
array were utilized for this purpose. This then provided a template for 


three subsequent copies of the Sum_amplitude process each of which was 
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individually altered to process a segment of the range indices, and then 
stored on a disk file. 

Similarly, the addresses of variables that were local to the 
Indata procedure were amended so they were within an SBC's private 
memory. Having loaded the Trace matrix from the floppy disk data file 
into the appropriate memory locations that had been assigned by the 
compiler, the entire Indata version (which included the data structures) 
was saved on a disk file. A subsequent reloading of this file into 
memory caused the Trace matrix to also be available and so a separate 
loading of the Trace matrix data file was not necessary. For the 
experimental demonstration of the feasibility of the multiple micropro- 
cessor concept the Indata process controlled the bearing parameter and 
caused the entire algorithm to scan ten complete sectors and then stop. 
An operational version would continue indefinitely. As the double buf- 
fering design resulted in two bearing increments being processed on each 
pass through the procedure, it was necessary to ensure that each sector 
consisted of an even number of bearing increments, i.e. thirty vice 
thirty-one. No attempt was made to include pulse propagation time or 
data acquisition time. The procedure that would fulfill these functions 
was left as a "stub" since it was not required for the demonstration of 
the synchronization between the processes. 

3. Results 

The capabilities of the MDS and DDT were used to subsequently 
load the Indata process into the private memory of one of the SBC's, and 
the Sum_amplitude processes into the private memory of each of the three 
remaining SBCs. Figure 18 provides a block diagram of the experimental 
configuration where process D is the Indata process and processes 24 to 
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Figure 18. Block diagram of experimental 
configuration. 
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S3 are the Sum amplitude processes. With the processes executing simul- 
taneously the time required for the amplitude summation function to be 
per formed on ten sector scans was recorded. 

Executing the algorithm with a different number of Sum_ampli- 
tude processes in parallel, while ensuring that all range indices were 
processed, permitted an evaluation of the effect of increasing the number 
of parallel processors. The resulting times for one sector scan and 
consequently the performance in terms of scans per second are given in 


Table 2 and Figure 19. 


Time for one scan Scans per 
Processes (seconds +0.01) Second 
BA Sy 6.45 0.155 
D, Sy> S9 Sas 0.310 
D, Sy: Sos S3 2.15 0.465 


TABLE 2. Results of the experimental demonstration of the 
feasibility of implementing the signal processing 
algorithm on multiple microprocessors. 

The correct operation of the algorithm together with the per- 
formance times indicated that the design was a feasible approach to 
obtaining near real-time execution of the proposed signal processing 
algorithm. The fact that the processes were stored in the SBCs private 
memory and required access to the bus only when the common data struc- 
tures and semaphores were referenced resulted in a minimum of bus usage. 
An examination of the code for the innermost loop of the Sum_amplitude 


process (the summation of the Indata array elements identified by the 
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trace) indicated that the bus usage was 7.9% of the loop's execution 
time. The number of parallel SBCs executing simultaneously could theo- 
retically be increased to twelve before bus contention would adversely 
influence the overall execution time. Table 2 and Figure 19 depicting 
the results of the experiment clearly indicate the linear relationship 
between performance, in terms of scans per second, and the increasing 
number of parallel processors used. If the number of SBCs operating in 
parallel were increased to seven, plus one for the Indata process for 

a total of eight, then a scan rate of one sector per second would be 


achieved. 
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¥. CONCLUSIONS 


For the problem of locating and classifying objects buried in 
marine sediments an underwater acoustic imaging system has been pro- 
posed and the thrust of this thesis has been the design of the signal 
processing algorithm that would extract the information from scattered 
acoustic signals for subsequent operator display. Using a developmen- 
tal form of the algorithm it was determined that the acoustic imaging 
system should be a narrow beam system for optimum resolution and 
imaging performance. Since the processing time was excessive, the 
algorithm was redesigned for implementation on multiple processors for 
execution efficiency. It was experimentally demonstrated that this 
design was a feasible approach to obtaining near real-time execution of 
the algorithm. 

The maximum rate of sector scanning that is possible for the pro- 
posed system is slightly in excess of two sector scans per second. 
This was arrived at by considering the pulse propagation time and the 
data acquisition time, a total of 14.5 milliseconds at each bearing, 
as the minimum time between pulses. It has been shown that the use of 
multiple microprocessors can result in an efficient implementation of 
the proposed signal processing algorithm. The number of parallel pro- 
cessors required to reduce the amplitude summation process to about 15 
milliseconds (and thus achieve the maximum rate of two scans per 


second) would be dependent on the microprocessor chosen. 
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VI. RECOMMENDATIONS 


The signal processing algorithm and its implementation on mul- 
tiple microprocessors, as discussed in this thesis, would form the core 
of a special purpose processor for the acoustic imaging system. Two 
areas would require development and detailed definition before the 
entire special purpose processor could be designed: data acquisition 
to provide inputs to the processing algorithm, and output from the 
algorithm in the form of a display. 

It is recommended that the Indata process outlined in the thesis 
be expanded from the current "stub" status to perform the functions 
previously mentioned; i.e. control of pulse transmission timing, 
sampling the hydrophone outputs, and storing the results into the two 
buffer arrays. It is anticipated that the approaches used for the 
algorithm design would be readily applicable to these functions. How- 
ever, the actual details would be highly hardware dependent which would 
require investigation prior to designing a parallel/pipeline approach to 
this portion of the specialized processor. 

Similarly, the interface between the multiple microprocessor 
system (performing the signal processing) and a display processor would 
require definition and would be hardware dependent. It is considered 
likely that the image information after each sector scan could be 
stored thus forming a three dimensional representation in memory (the 


third dimension resulting from the travel of the underwater vehicle 


os 





over a section of sediment). This would provide an opportunity for a 

two dimensional display in either the vertical plane (such as range and 
bearing) or in a horizontal plane at some range. Furthermore, the possi- 
bility of a three dimensional display could be considered. It is evi- 
dent that there is much promising work to be accomplished in order to 


produce the special! purpose processor for the acoustic imaging system. 
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